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We consider the effect of a positive cosmological constant on spherical gravitational collapse to a 
black hole for a few simple, analytic cases. We construct the complete Oppenheimer-Snyder-deSitter 
(OSdS) spacetime, the generalization of the Oppenheimer-Snyder solution for collapse from rest of 
a homogeneous dust ball in an exterior vacuum. In OSdS collapse, the cosmological constant may 
affect the onset of collapse and decelerate the implosion initially, but it plays a diminishing role as 
the collapse proceeds. We also construct spacetimes in which a collapsing dust ball can bounce, or 
hover in unstable equilibrium, due to the repulsive force of the cosmological constant. We explore 
the causal structure of the different spacetimes and identify any cosmological and black hole event 
horizons which may be present. 
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I. INTRODUCTION 

Recent measurements of Type la supernovae suggest 
that our universe may have a nonzero cosmological con- 
stant A > |I],p| . A more recent analysis || of the pecu- 
liar motion of low-redshift galaxies seems to give further 
evidence for a finite A. Obviously, such an interpretation 
of the data, if correct, will have huge implications for cos- 
mology. More generally, if a cosmological constant must 
be restored to Einstein's equations of general relativity, 
surprises may turn up in other physical applications of 
Einstein's field equations, although the small size of the 
constant precludes its having a significant effect on the 
scale of typical galaxies, stars or planets. It is therefore 
interesting to consider, at least as a point of principle, 
what impact, if any, the presence of a finite cosmological 
constant has on our conventional picture of gravitational 
collapse to a black hole. Many of the important dynami- 
cal and geometric features of catastrophic collapse in the 
absence of a cosmological constant are revealed by the 
analytic Oppenheimer-Snyder model Q , which describes 
the collapse of a spherical, homogeneous dust ball, ini- 
tially at rest in an exterior vacuum, to a Schwarzschild 
black hole. In this paper we generalize this model ac- 
counting for the presence of a positive cosmological con- 
stant. We also consider closely related, dust ball solutions 
for which the implosion does not begin at rest. 

There are a number of questions that motivate our 
analysis: How does the cosmological constant, which acts 
as a repulsive force, affect the motion and fate of a col- 
lapsing object? Under what circumstances, if any, can a 
cosmological constant prevent the collapse of a dust ball 
which is initially imploding? What is the global horizon 
structure of an exponentially expanding universe contain- 
ing a collapsing dust ball ? When do black holes form? 



II. DYNAMICS OF A HOMOGENEOUS DUST 
SPHERE WITH COSMOLOGICAL CONSTANT 

The interior of a homogeneous sphere is given by the 
Friedmann-Robertson- Walker (FRW) metric 



ds 2 = -cIt 2 



dx 2 



1 — kx 2 



(2.1) 



where k — —1, or 1 for a hyperbolic, flat or spheri- 
cal spatial geometry, respectively. The density remains 
homogeneous on spatial slices of constant time r. The 
surface of the sphere is located at constant x — X, where 
< X < 1. Einstein's equations for a pressureless fluid 
(i.e., "dust") of density fi/a 3 in the presence of cosmo- 
logical constant A yield 
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This interior spacetime is often called the Friedmann- 
Lemaitre (FL) universe. 

The standard Oppenheimer-Snyder (OS) solution Q] 
for the interior of a collapsing homogeneous dust sphere, 
initially at rest, is a piece of a closed FL (k — 1) uni- 
verse with A = 0. In OS collapse, the initial time-slice 
t = is defined at the moment of time-symmetry at max- 
imum expansion, when the right hand side of Eq. (2.2) 
vanishes. The same form for the interior metric with 
k = 1 applies in the presence of a cosmological constant. 
The solution we seek - collapse of a spherical dust ball 
from rest in an exponentially expanding universe with a 
positive cosmological constant - we shall refer to as an 
Oppenheimer-Snyder-de Sitter (OSdS) spacetime. 

According to the generalized Birkhoff theorem JS) , the 
vacuum spacetime outside the sphere is Schwarzschild-dc 
Sitter (SdS) § 
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where M is a constant. This spacetime represents the 
vacuum exterior of a spherical mass immersed in the ex- 
ponentially expanding (at large r) de Sitter space. 

We neglect for the moment t he presence of the dust 
sphere and assume that metric (2.3) describes the entire 
spacetime with r in the range < r < oo. The metric 
function f(r) (see Fig. 1) reaches the maximum value 

/ max = 1 - (9M 2 A) 1/3 at r m = (3Af/A) 1/3 . Thus, for 
M < 1/3A 1 / 2 , f(r) has two real positive roots, rh and 
r c > rh, where 
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(we choose < sin" 1 A < it/2 for < A < 1; the third 
root is negative, r3 = — r^ — r c ). For all null and timelike 
geodesies crossing the black hole horizon r = inward, 
— d r is the future-directed timelike vector, and so they 
all terminate at the singularity r = 0. All the future 
and outward-directed (timelike or null) geodesies at r > 
rh cross the cosmological horizon r = r c and ultimately 
reach cosmological null infinity I + at late times. 

Returning now to the motion of a dust sphere with the 
exterior spacetime (2.3), we match the spacetime geome- 
try across the sphere's surface r = R(t) (x = X = const) 
which requires (1) the continuity of the surface's 3-metric 



(3) ds 2 = -dr 2 + R 2 (dO 2 + sin 2 



(2.5) 



where R(t) = Xa(r), as well as (2) the junction condi- 
tion \K)] = K]\ t —K)\. — 0, for the extrinsic curvature 



K, 
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(2.6) 



The three vectors are intrinsic to the spacetime hy- 
persurface swept by the moving surface, u = tdt + Rd r 
(t = dt/dr), dg and d^, while n = n*<9 t + n r d r is the 
outward-directed unit 4- vector orthogonal to the surface, 
erom the orthonormality conditions 
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and the metric forms ( |2.l| ) and (^_^), we obtain the 
t— component of the 4- velocity 
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FIG. 1. The effective potential f(r) (A = 0.07/M 2 , 
X — 0.96). Oppenheimer-Snyder-de Sitter (OSdS) collapse 
starts from rest at the dot and proceeds inward (leftward 
along the dotted line). 



and components of the unit normal 

/ . \ 1/2 

n r = m(f + R 2 ) , 

(r)2 = ±1) on the outside, and 
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n T = 0, 
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on the inside of the surface. 
The condition [Kg] = 0, where 
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implies 7/2 = 1 and leads to the equation that describes 
motion in the effective potential f(R) (see Fig. 1) 



R + f(R) = l-kX 2 



Comparing Eqs. (|2.2|) and (2.12) we identify 



M = ^„A 3 
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(2.12) 
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(The junction conditions for other components of the ex- 
trinsic curvature do not yield additional information.) 
For k = 1 there is a range of parameters, M < 1/3A 1//2 

and X > (9Af 2 A) 1/6 , for which equation = 1-A 2 -/ = 
Ar 2 /3 + 2A//r — A 2 has three real roots, out of which two, 
r and rb > r are positive 
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FIG. 2. A dust ball with an exterior whose spac elike hyper- 
surfaces have cylindrical geometry [see Eq. (2.15)]. One an- 
gular dimension is suppressed. The dust ball recedes leftward 
(y < 0) as viewed from the exterior. The interface (x — X) 
normal n is drawn tangential to the sphere's interior. 



The third root, r% = — — r a is then negative. Accord- 
ing to Eq. (2.12), the sphere at R = r Q or R = r\, is 



momentarily at rest. 

For f{r) > t is a timelike future-directed coor- 
dinate and Eq. J2.8| ) implies 771 = 1. If, on the other 
hand, r < or r > r c , so that f(r) < [assuming f(r) 
has positive roots], r is timelike. If, in addition, — d r is 
future-directed [61, it is convenient to introduce the new 
notation, T = —r < and y = t, in which metric ( |2.3| ) 
takes on the time-dependent form 



ds 2 = --dT 2 +gdy 2 + T 
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Metric (2.15) describes a 3-dimensional cylindrical space- 
like hypersurface of radius \T\ , symmetric around its y- 
axis (see Fig. 2). For r < rt, this hypersurface contracts 
radially and expands in the y-direction from zero size at 
the coordinate singularity T = — r^, to an infinite extent 
at the physical singularity T = 0. For r c < r < 00, it 
contracts in the radial direction and along its y-axis un- 
til the coordinate singularity g(T) = is reached at time 
T = -r c . 

The vacuum spacetime given by metric ( |2.15| ) may ex- 
tend to infinity in both directions along the j/-axis. Alter- 
natively, it can be bounded on the left (so that n points 



toward increasing y, 



> 0; see Fig. 2) by the sur- 



face of an inevitably contracting sphere [R = — T < 0, 
and thus 771 = —1; see Eq. (2.9)] that recedes leftward, 
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along the y axis. The trajectory of the surface is given 
by 
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For r < Hi, the sphere's surface therefore starts at T — 
— rh, where g oc T — — * 0, at which point y — * —00 
and y — > 00, and reaches some finite y = yf at T = 0. 
For r > r c , the surface travels from a finite y = yi, at 
T = —00, to y — — 00 at T = —r c . 

If, on the other hand, we place the contracting spher e 
to the right of the vacuum sector given by metric ( 2.15 ) , 
n y < (and thus 771 = 1), and the sphere's surface will 
recede toward increasing y, y > 0. If r is a future- 
directed timelike coordinate |l(j (so that the notation 
change would be T = r, y = t), the motion of the sphere 
is reversed: it expands, R = T > and, if, e.g., placed 
on the left, n v > 0, it advances to the right, y > 0. 



III. BOUNCING SPHERE 

As discussed in the previous section, for k = 1 we 
have the range of parameters, M < 1/3A 1 / 2 and X > 

(9M 2 A) 1 ^ 6 , for which a dust sphere contracting from 
large radii will bounce at R = > t\, where 
is the larger of the two positive roots of the equation 
f(r) = 1-X 2 (see Fig. 1). 

The Penrose diagram for the complete spacetime is 
given in Fig. 3. At large negative time r (between the 
past null infinity I~ and the past cosmological hori- 
zon Tl~) the dynamics is dominated by the cosmologi- 
cal constant and the sphere contracts exponentially, R oc 
exp(— y/AjEr). After the sphere's surface has crossed 

it will reverse its motion at R = Xa{r = 0) = 
and then expand toward the future cosmological hori- 
zon and ultimately the future null infinity I + . In 
contrast to an implosion in a static Schwarzschild back- 
ground, the presence of a positive cosmological constant 
is sufficient to halt and reverse the collapse in this case. 
Here the presence of the dust does not alter in essence 
the familiar bounce of the spherical spatial hypersurfaces 
near the "throat" of the de Sitter spacetime [ 111 . 

The presence of a massive sphere, however, is not en- 
tirely without an effect on the global structure of the 
spacetime. In the limit of vanishing mass M, the ra- 
dial null ray that enters the sphere at the exact moment 
when its surface crosses the future cosmological horizon 
(7i+)radius, R = r c , will take an infinite time r to reach 
the re-expanding sphere's center x — 0. On the other 
hand (see the Appendix), for any finite mass M > 0, the 
null ray will reach the center in a finite time, then pro- 
ceed outward and reach some finite x = Xf at infinitely 
late times. In this case, the surface r c no longer has the 
significance of a cosmological horizon |T^] with respect to 
an observer at x = 0. 
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FIG. 3. The Penrose diagram for the spacetime containing 
a "bouncing" dust sphere. The bouncing dust sphere to the 
right (identical to the one on the left side) is introduced for 
geodesic completeness. 

The spacetime shown in Fig. 3 is geodesically complete, 
i.e., all geodesies extend to infinite values of their affinc 
parameters in both directions. To achieve the complete- 
ness we have placed another dust sphere at the right of 
Fig. 2 and Fig. 3. While the second sphere must have 
the same values of M and A, we are free to choose for 
it any values of k r and X T , independent of the choices k 
and X for the sphere at the left. Of the great variety of 
combinations of two dust spheres, we will in this paper 
discuss only a few. The choice we have made in Fig. 3 is 
symmetric, i.e, fc r = k = 1 and A r = X. Since the cen- 
ters of the spheres cross the null rays 7i+ or TL~ , the two 
spheres are in causal contact. [Other ways of extending 
the spacetime to the right — with either massive spheres 
or interiors of black holes — are discussed in the next 
section.] 

IV. COLLAPSING SPHERE 

For the range of parameters discussed in the previous 
section, M < 1/3A 1 / 2 and X > (9Af 2 A) 1/6 , one can 
follow the motion of the sphere at R < Ri, where Ri is 
smaller of the two positive roots of the equation /(r) = 
1 — X 2 . The sphere's surface springs out of the past 
singularity at r = (see Fig. 4), then emerges through 
the past black hole horizon TC^ , reverses its expansion at 
R = Ri, plunges through the future black hole horizon 
Wfc", and finally ends in the future singularity at r = 0. 
As discussed in Section 2, the spacetime exterior to the 
sphere, usually called the Schwarzschild-de Sitter (SdS) 
spacetime ||, is characterized by two pairs of horizons: 
the future and past black hole horizons at r = r^, and the 
future Tic and past 7i~ cosmological horizon at r = r c . 

In this scenario the sphere always stays inside the cos- 
mological horizon at r c . If we follow the sphere's implo- 
sion from the moment of time-symmetry at maximum 



FIG. 4. Collapse to a black hole in the Schwarzschild-de 
Sitter universe. The dust sphere's surface emerges from the 
past singularity through the past black-hole horizon, stops 
its expansion at point ^ and then recollapses to a black hole. 
The Oppenheimer-Snyder-de Sitter (OSdS) collapse to a black 
hole starts from the moment of time-symmetry. The bouncing 
sphere, X r > (9Af 2 A) 1 ^ 2 , of mass M on the right side is in- 
troduced for geodesic completeness; it can be replaced by an- 
other re-collapsing sphere. The black-hole horizon structure 
corresponds to region II in the parameter space of Fig. 6. 

expansion R = r Q in Fig. 4), the motion is a straight- 
forward generalization for a non-vanishing cosmological 
constant A of the familiar Oppenheimer-Snyder collapse. 
The collapse is homologous and the density remains ho- 
mogeneous on t = constant time slices. Qualitatively, 
the cosmological constant serves as a perturbation whose 
influence on the collapse diminishes as the collapse pro- 
gresses. We illustrate the initial braking effect of A on 
the collapse in Fig. 5 for the indicated values of AM 2 . 
We measure the sphere's time in the units of the proper 

1 /2 

time, T\ =0 = it (i2f/8MJ , that a sphere would take 
to evolve from the static point at R = Ri all the way to 
the final singularity in the case of a vanishing A. 

In the left portion of Fig. 4, the past black hole horizon 
TC^ reaches the center of the sphere x = after the fu- 
ture black hole horizon 7ijJ~ has emerged from the origin. 
This configuration of the black-hole horizons would, e.g., 
make it impossible for an observer at x = to take off in 
the radial direction and escape from the black hole after 
receiving the earliest possible signal from the past null 
infinity I~ . As we discuss in the Appendix, this is the 
case for the sector II in the plane of parameters X and 
(9AM 2 ) 1/6 shown in Fig. 6. 

For the values of the parameters from sector I of Fig. 6, 
on the other hand, the arrival of TC^ at x — precedes 
the departure of Ht. This case is shown in Fig. 7. 

For any values of M and A (including A = 0) that allow 
a recollapse (or the OSdS implosion from rest), M < 
1/3A 1 / 2 , the earliest null ray emitted from the portion 
of the past naked singularity not covered by dust (see 
Figs. 4 and 7) will reach the sphere's center before the 
emergence of Hjj" (see the Appendix) . This will allow an 
alert observer at x = to avoid the future black-hole 
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FIG. 5. OSdS collapse to a black hole. Shown is the evolu- 
tion R(t) of the surface of a dust sphere starting from rest at 
R = 3M for the indicated values of AM 2 . The solid dots mark 
the points at which the spheres' surfaces cross their black-hole 
horizons. 



FIG. 7. A recollapsing sphere (compare Fig. 4) for param- 
eters from region I of Fig. 6. Instead of another sphere to 
the right, the spacetime is closed off with an infinite series of 
alternating Schwarzschild-de Sitter exteriors and black-hole 
interiors. 




FIG. 6. Regions in the parameter space that allow an ob- 
server to escape from the recollapsing sphere to infinity after 
receiving at x — the first signal from the dust-free portion 
of the past naked singularity (I+II) or from the past null 
infinity (I). 
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FIG. 8. A static dust sphere in unstable equilibrium for 
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singularity if he heeds the warning coming from visible 
(assuming light can get through any finite-density region) 
sector of the initial singularity. 

Instead of another sphere shown in Figs. 3 and 4, Fig. 7 
is completed at the right by an infinite sequence of alter- 
nating Schwarzschild-de Sitter exteriors (with their char- 
acteristic cosmological horizons and spacelike null infini- 
ties) and black-hole interiors (bounded in the past and 
future by black-hole singularities). 

For the special value X — (9M 2 A)^ 6 , the sphere can 
hover in a state of unstable equilibrium at the maximum 
of the effective potential at R — r m . The Penrose dia- 
gram containing such a sphere is shown in Fig. 8. 

Finally, for M < 1/3A 1 / 2 but X < (9M 2 A) 1/6 , a 
sphere contracting from i? = oo at r = — oo will pass 
through 7i~ and then form a black hole. Fig. 9 contains 
the corresponding Penrose diagram. Notice that at early 
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FIG. 9. Collapse into a black hole for X < (9Af 2 A) 1/6 . 
The dust sphere's surface contracts through the past "cosmo- 
logical horizon" He , and then collapses to a black hole. 



times the spacetime structure is similar to the early-time 
portion of the bouncing case (see Fig. 3), while the late- 
time evolution is akin to that of Fig. 4. 



V. GLOBAL COLLAPSE 

If M > 1/3A 1 / 2 , f(r) is negative everywhere and there 
is no static (i.e., with a time- like Killing vector) portion of 
the spacetime. The entire spacetime is analogous to the 
familiar vacuum Schwarzschild solution inside the event 
horizon, where the spacetime is dynamic and the radial 
and time coordinates reverse roles [see the last three para- 
graphs of Section 2] . Metric ( 2.15| ) describes a cylindrical 
spacelike hypersurface that shrinks (if the timelike coor- 
dinate T = — r is future directed) in the radial (T > 0) 
direction and expands along its y-axis {y = t) as the 
singularity at T = is approached. At the left edge of 
this vacuum spacetime (see Fig. 10), the surface of the 
collapsing spher e rec edes leftward along the y-axis, ac- 
cording to Eq. ( ^.16| ). 

Since the integral J dy = — J a oQ dT/g(T) is finite, a 



null ray travels only a finite difference Ay [see Eq. (2.17)] 
in the spacetime outside the sphere over the entire evolu- 
tion from the past null infinity, T = — oo, to the singular- 
ity at T = 0. This allows us to to place another identical 
collapsing sphere (needed for geodesic completeness) re- 
ceding to the right (y > 0) of the first one (see Fig. 10) 
so that the two spheres either have causal contact (the 
specific case of Fig. 10) or are causally disconnected. In 
either case, the entire spacetime ends in a global cosmo- 
logical singularity. 

For k = or k = — 1, the static, momentarily static 
and bouncing solutions described above do not exist. De- 
pending on whether M < 1/3A 1 / 2 or not we have, 
respectively, a black hole formation shown in Fig. 9 or a 
"big crunch" shown in Fig. 10. 



VI. SUMMARY 
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FIG. 10. Collapse into a global (cosmological) singularity 
for M > 1/3A 1/2 . 



and rebound en route to an exponential expansion at 
late times. In contrast to the previous case, the bounc- 
ing sphere's evolution is dominated by the cosmological 
constant, with the matter playing the role of perturba- 
tion. 

For all values of the parameters outside the above 
range, an initially contracting sphere will not be able 
to avoid being infinitely squeezed in the final singularity. 
The fate of the exterior spacetime depends, however, on 
the mass of the sphere: for M < 1/3A 1 / 2 , the sphere 
forms its own black hole within its horizon, allowing the 
exterior space to expand exponentially at late times. For 
M > 1/3A 1 / 2 , the sphere drags the entire spacetime into 
a "big crunch": the exterior contracts on its way to the 
familiar de Sitter-like "throat" (see, e.g., pj]] ) but can- 
not escape out of it due to the overwhelming pull of the 
sphere's gravity. 
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In this paper we have investigated the influence of a 
finite cosmological constant A on the evolution of a ho- 
mogeneous sphere made of pressureless matter. In ad- 
dition to A, the evolution is determined by the mass M 
of the sphere and the comoving extent, X, of the sphere 
interior's spherical (fc = 1), flat (k = 0) or hyperbolic 
(k = —1) slices of homogeneity. 

A straightforward generalization of the familiar 
Oppenhcimer-Snyder collapse from rest is possible only 

for M < 1/3 A 1 / 2 , k = 1 and X > (9M 2 A) 1/6 . In this 
case, the cosmological constant can slow down the col- 
lapse initially, but at later times the sphere's self-gravity 
dominates entirely and eventually pulls the sphere into 
the final singularity. 

The same range of parameters allows, however, a 
sphere contracting from large radii to avoid a collapse 



APPENDIX A: EXTENSION OF THE 
COSMOLOGICAL AND BLACK HOLE 
HORIZONS INSIDE THE DUST SPHERE 

1. Bouncing sphere 



Solving Eq. (2.1) for future- and inward-directed ra- 
dial null rays along 7i+ and integrating from the point 
(R = r c ) at which the re-expanding dust sphere's surface 
crosses 7i+ (see Fig. 2), we obtain 
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F(X,p), 



(Al) 



where p = Ar 2 /3 = 1 — 2M/r c , z = R/r c and X{ is 
the outermost point inside the dust sphere that can be 
reached by the null ray after it started at x — X and then 
passed through the center x = 0. In Eq. (Al) we have 
introduced the conformal time £(r) of the FRW metric. 

In the limit of vanishing mass, p(M = 0) = 1, 
F(X, 1) = sin -1 X and thus Xf = 0. This is the famil- 
iar structure of the de Sitter spacetime |l2). For M > 0, 
however, Xf > 0, and thus the center of the sphere, x = 0, 
will cross 7i+ at a finite proper time t c (as shown in 
Fig. 2). 

This can be shown as follows. For any given A, 
(dr c /dM)\ = 2/(1 — 3p). Since p has the minimum 
value p = 1/3 for M = 1/3A 1 / 2 (r c = 3M = 1/A 1 / 2 , 
/ max = 0), (dr c /dM) A < and thus (dp/dM) A will 
always be negative. It then follows (dF / dM)x = 
(dF/dp) x (dp/dM) A > 0. Hence, for any M > 0, 
F[X,p(M)] > F[X,p(M = 0)] = sin" 1 A and thus 
Xf > 0. The surface r — r c , therefore, loses for M > 
the significance of an event horizon (with respect to the 
world line x = 0) that it has at M = 0. 



C+ = 2C(r ) - C(r h ) 
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If the past horizon is to reach the sphere's center before 
the emergence of the future horizon, we must therefore 
have 



C(r„) > C(nO 



sin -1 A, 



(A5) 



the condition satisfied by all points in region I of the pa- 
rameter space of Fig. 5 and by the specific case shown in 
the Penrose diagram of Fig. 6. If, on the other hand, we 
require that the null ray, emanating from the boundary 
between the sphere and the dust-free portion of the past 
singularity, gets at C = sin -1 A to the sphere's center be- 
fore the emergence of Ti.^ , the following inequality needs 
to hold 



1 



C(r„) ~ hCK) > sin" 1 A. 



(A6) 



It turns out that condition (A6) is satisfied for all values 
of the parameters that allow a recollapse or an OSdS 
implosion. 

In the case of a vanishing cosmological constant, 



C( r o) ~ C( r h) = 7T — 2 sin A and the condition (A5) is 



satisfied if A < v3/2- At the same time, C( r o) — ^C( r ii) 
7r — sin -1 A, and the inequality ( |A6| ) holds if A < 1, 
which is satisfied by all recollapsing spheres. 



2. Collapsing sphere 



We now turn to the configuration of the past (H h ) 
and the future (7~U) black hole horizons. Measuring the 
conformal time of the FRW metric from the initial sin- 
gularity 



R ' x da 



dR 



o \R(r - R)(r h - R)(fa 



i?)| 1/2 ' 
(A2) 



where A = 9AM 2 , and again solving Eq. (2.1) along the 
past black-hole event horizon, we find the value 



C- = C(nO + 



sin" 1 A, 



(A3) 



of the conformal time when the past black-hole horizon 
reaches the center of the sphere (see fig. 3). Notice 
that the conformal time depends only on A and AM 2 in 
addition to the current value of R. 

Since the sphere's evolution is time-symmetric around 
the point of reversal r , the conformal time that elapses 
between the initial and past singularities is 2£(r ). 
Hence, the future black hole horizon 7ijJ" emerges fron 
x = at the conformal time 
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